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Resume 

La formalisation non standard, due a P. Cartier et Y. Perrin, des oscillations rapides 
fournit un cadre mathematique adequat pour de nouvelles techniques d'estimation 
non asymptotiques, ne necessitant pas l'analyse statistique habituelle des bruits 
entachant tout capteur. On en tire diverses consequences sur les bruits multiplicatifs, 
la largeur des fenetres d 'estimations parametriques et les erreurs en rafales. Pour 
citer cet article : M. Fliess, C. R. Acad. Sci. Paris, Ser. I 342 (2006). 

Abstract 

Noises: a nonstandard analysis. Thanks to the nonstandard formalization of fast 
oscillating functions, due to P. Cartier and Y. Perrin, an appropriate mathemati- 
cal framework is derived for new non- asymptotic estimation techniques, which do 
not necessitate any statistical analysis of the noises corrupting any sensor. Various 
applications are deduced for multiplicative noises, for the length of the parametric 
estimation windows, and for burst errors. To cite this article: M. Fliess, C. R. Acad. 
Sci. Paris, Ser. I 342 (2006). 
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Abridged English version 
1 Introduction 

Recent works (see, e.g., [5,6,7,8,9,10,11,12]) in signal processing and automatic 
control yield quite efficient estimation techniques, which are non-asymptotic 
and do not necessitate any statistical treatment of the noises corrupting any 
sensor. An appropriate mathematical framework is provided by the nonstan- 
dard formalization, which is due to Cartier and Perrin [2], of fast oscillating, 
or fluctuating, functions. Various applications are derived for multiplicative 
noises, for the lengths of the estimation windows, and for burst errors. 



2 Noises 

Let *N, *M be the usual nonstandard extensions of N, M. Replace [0, 1] 1+D C 
R 1+D by the hyperfinite set I = I x L x . . . x I D , I t = {0, ^, . . . , &=±, 1}, 
where N t G *N, i = 0,1,..., D, D G N, is unlimited. The' function \ : 
I t \{l} — > *R, j±- h-> < k L < N L - 1, is the Lebesgue measure of 1 L [2]. 
The Lebesgue measure A = Ao x Ai x . . . x of I is clear. See [2] for the 
notions of quadrable sets and of Lebesgue integrability. A S'-integrable function 
/ : I — > *R is fast oscillating [2], or fast fluctuating, if, and only if, any integral 
f A fal\, where A C I is quadrable, is infinitesimal. Let m : I — > *R be a 
Lebesgue integrable function. A noise, of mean m, is a S'-integrable function 
n : I — > *R such that, Ml — 0,1, . . . ,D, M ^ k G I k , k 7^ the projection 
n (^o, Ci, • • • ) £d) — "i(Co) Ci> • • • i Cr>) : I t — >■ *K is fast oscillating. If it is possible 
to take m = 0, n is said to be zero-mean. 



3 Applications 

3.1 Multiplicative and additive noises 

Let x : I — > *M be a function of class S l [3], which we call a signal. Let the 
sensor y = ttix + n 2 be corrupted by a multiplicative noise tii, of mean 1, and 
an additive noise ri2, of arbitrary mean. 

Proposition 3.1 y = x + n where n = n 2 + (tii — l)x is a noise of the same 
mean as ri2. 
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3.2 Parametric estimations 



Set D — 0. Assume that the the shadow of the signal x : I — > ffi is analytic. 
Assume also that the constant parameter e M is linearly identifiable [6,10]. 
It leads to the estimator (4) where [0] e (£) is the estimate of 9 at £; [0, t] is the 
estimation window of length t; the shadow of is an analytic function 
[0, 1] — > M, called divisor, such that 5(0) = 0; n is an additive noise. 

Proposition 3.2 Assume that n is zero-mean. If the length of the estimation 
window does not belong to the monads of the divisor's zeros, the estimate 
[0] e (t) belongs to the monad of 9. The length is then appreciable. If the length 
is infinitesimal, [6] e (t) does not belong necessarily to the monad of 9. 

The parameter, often called symbol, which has to be demodulated (see, e.g., 
[17,21]) is generally associated to a signal satisfying a linear differential equa- 
tion with time-polynomial coefficients. Burst errors may be understood in our 
setting in the following way: contrarily to Proposition 3.2 the noise is not 
zero-mean. Since the transmission duration of any signal is short, we may as- 
sume that the shadow of the unknown mean m(t) is a polynomial p(t) of a 
given limited degree. Consider p as a structured perturbation which may be 
annihilated [6,10] by some given limited power of ^. 

Proposition 3.3 Assume that the shadow of the signal is not annihilated by 
a limited power of There exists then an estimate [9] e (t) of the symbol which 
belongs to the monad of 9, for any t which does not belong to the monad of a 
divisor's zero. 



1 Introduction 

Des travaux recents (voir, par exemple, [5,6,7,8,9,10,11,12]) conduisent a des 
techniques d'estimation non asymptotiques, efficaces en signal et en automa- 
tique, sans aucun recours aux traitements statistiques habituels des bruits 
entachant tout capteur. Rappelons brievement de quoi il retourne pour ces 
methodes, illustrees par maints exemples, mal maitrises par ailleurs. On dis- 
tingue [6,10,11,12] deux types de perturbations, celles dites structurees, anni- 
hilees avec des operateurs differentiels lineaires, et celles dites non structurees, 
considerees comme des oscillations, ou fluctuations, rapides, attenuees par des 
nitres passe-bas, comme les integrates iterees. 

La formalisation non standard [23] des oscillations rapides, due a Cartier et 
Perrin [2], generalisant des travaux anterieurs de Harthong [14] et Reder [22], 
fournit un cadre mathematique, permettant d'introduire les bruits et leurs 
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moyennes dans un cadre entierement deterministe. Plusieurs justifications is- 
sues de la pratique des ingenieurs sont proposees. On en deduit divers resultats 
sur les bruits multiplicatifs, la largeur des fenetres d'estimations, que confor- 
tent deja de nombreuses simulations numeriques (voir [5,6,7,8,9,10,11,12] et 
leurs bibliographies) et quelques experiences de laboratoire, et sur les erreurs 
en rafales, burst errors en americain, pour lesquelles ce travail reste a faire. 



2 Bruits 

Nous avons redige selon le formalisme ZFC de Robinson [23] , tout en emprun- 
tant beaucoup au vocabulaire de [2] et [3], qui se situent dans VISTde Nelson 
[18]. 



2.1 Definition non standard 

Designons par *N, *R les extensions non standard de N, R. Remplagons [0, 1] 1+D C 
R 1+D par l'ensemble hyperfini I = I x Ii x . . . x I D , l b = {0, ±-, . . . , Zj=±, 1}, 
ou N L G *N, i = 0, 1, . . . , D, D G N, est illimite. La fonction A t : I t \{l} U *R, 
j±- i— > jjr, < k L < N t — 1, est appelee [2] mesure de Lebesgue de I t . On 
definit alors, de fagon evidente, la mesure de Lebesgue A = AoxAiX...xAd 
de I. Renvoyons a [2] pour les notions d'ensembles quadrables et de Lebesgue- 
integrabilite. Une fonction S'-integrable / : I — > *K. est dite a oscillations, 
ou fluctuations, rapides [2] si, et seulement si, toute integrate J A fdX est in- 
finitesimale, ou A C I est quadrable. 

La definition suivante paraphrase le theoreme 9.3.6 de [2] sur la decomposition 
d'une fonction .S-integrable en somme d'une fonction Lebesgue-integrable et 
d'une fonction S'-integrable a oscillations rapides. Soit m : I — > *R une fonction 
Lebesgue-integrable. Un bruit, de moyenne m, ou, en adaptant la terminolo- 
gie expressive de [3], crepitant autour de m, est une fonction S'-integrable 
n : I — > *R telle que, V t = 0,1,..., D, V £ K G l K , k ^ i, la projection 
n (^o, Ci, • • • ) £d) — m(£o, Ci) • • • ) £d) '■ L *R est a oscillations rapides. La 
moyenne n'est, evidemment, pas unique : m' est aussi une moyenne, si, V i — 

0,l,...,D,V^eI K ,K^t, la projection m' (^o,6, • • • ,£i0-to(£o,6l, • • .,&) ■ 
L — > *R est infinitesimale presque partout (voir la proposition 9.3.12 de [2]). 
Si Ton peut prendre m = 0, on dit que n est de moyenne nulle, ou centre. Ce 
qui suit est aise : 

Proposition 2.1 Soit (ft : I — > *R fine fonction de classe S 1 [3j. Alors, le 
produit <ftn est encore un bruit, de moyenne <ftm. Si n est centre, (ftn Vest aussi. 
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2.2 Justifications 



2. 2. 1 Hautes frequences 

Soit n = Sfinio Asm.{VLt + </?), A, Q, </?, t e K. L'integrale /// n(r)dr, t i} tf £ K, 
est « petite » avec de « hautes » frequences f2. 

Remarque 1 Par contre, l'integrale Jtf (n(r)) 2 dT n' est pas « petite ». L'ecart- 
type d'un bruit, au sens du § 2.1, qui est, avec sa moyenne, de carre integrable, 
est done, en general, appreciable. 



2.2.2 Bruits blancs 

Considerons, avec bien des ouvrages pour ingenieurs en traitement du signal 
(cf. [21]), le capteur y d'un signal x bruite additivement 

y{t) = x{t) + n{t) (1) 



t G [0, 1], ou n(t) est un bruit blanc, souvent suppose stationnaire, centre et 
gaussien. Le lien avec le § 2.1 decoule de l'analogue echantillonne, usuel chez 
les praticiens (cf. [21]), 

y(aAt) = x(aAt) + n(aAt) (2) 



a = 0,l,...,iV, N e N, ou 

- At — jj est le pas d'echantillonnage, limite et appreciable, 

- les n(aAt) sont des variables aleatoires independantes, centrees, d'ecarts- 
types normalises a 1. 

En effet, toute somme 

*~ E n(aAt) (3) 

0<tj<aAt<t F <l 



tend presque surement vers avec le pas d'echantillonnage. 



Remarque 2 Les travaux plus mathematiques, comme [15], substituent a (1) 
V equation differentielle stochastique dy = xdt + dw, ou w est un processus 
de Wiener. La reecriture y(t) = y(0) + J * x(r)dr + w(t) — w(0) confirme, si 
besoin est, que Von ne peut se ramener a (1). Cette dichotomie sera examinee 
prochainement, de meme que ses liens avec certaines questions de mecanique 
stochastique (cf. [19,20]), en nous inspirant des moyennes glissantes de [14,22] 
et de la derivation quantique de [1]. 



5 



3 Applications 

3.1 Bruits additifs et multiplicatifs 

Soit x : I — > *R une fonction, appelee signal, supposee de classe S 1 , ayant, 
par consequent, une ombre C 1 [3]. Le capteur y : I — > *R de x est dit bruite 
additivement et multiplicativement si, et seulement si, y = tiix + ri2, oil les 
bruits multiplicatif ni et additif ti2 sont de moyennes respectives 1 et quel- 
conque, supposee, souvent, nulle. Le resultat suivant, qui repose sur le fait 
que (rii — l)x est, d'apres la proposition 2.1, un bruit centre, demontre que 
Ton peut se ramener au cas purement additif. On repond ainsi a des questions 
courantes de la litterature appliquee (c/. [13,16]) sur la maniere de traiter les 
bruits multiplicatifs. 

Proposition 3.1 On peut ecrire y = x + n, ou n = ri2 + (tii — l)x est un bruit 
de meme moyenne que n 2 . 

3.2 Estimations parametriques 

Nos estimations pour les images et les videos etant faites par balayage unidi- 
mensionnel [5], il est loisible de poser D = 0. Au vu des nombreux exemples de 
la litterature (cf. [21] et [5,6,7,10]), on suppose que le signal x possede une om- 
bre dans C w ([0, 1],M). Soit 9 G K. un parametre constant, suppose lineairement 
identifiable [6,10]. L'adaptation des calculs de [6,10] permet d'ecrire un esti- 
mateur sous la forme 



- t G Io, c G R, v > 0, k > 1 sont limites, 

- [0,t] est la fenetre d 'estimation, de largeurt, 

- 5(t) a pour ombre une fonction de C"^([0, 1], R), appelee diviseur, nulle en 



- n est un bruit additif, 

- [0} e (t) est l'estimee de 9 en t. 

3.2.1 Largeur des fenetres 

Supposons n centre. On deduit de (4) et de la proposition 2.1 le resultat 
suivant qui corrobore le caractere non satistique et non asymptotique, mais 





ou 
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non instantanee, de notre estimateur : 

Proposition 3.2 Si la largeur de la fenetre d 'estimation n'appartient pas au 
halo d'un zero du diviseur, I'estimee [6] e (t) appartient au halo de 9. Cette 
largeur est, alors, appreciable. Si, par contre, cette largeur est infinitesimale, 
[0] e (t) n'appartient pas necessairement au halo de 0. 

Remarque 3 L 'experience pratique « prouve » qu'une « courte » largeur suffit, 
en general, pour obtenir une bonne estimation. 

Remarque 4 Supposons que le bruit soit stationnaire et blanc. Sa variance 
etant proportionnelle a la distribution de Dirac a Vorigine, on dit (cf. [21]) 
qu'il contient toutes les frequences et que des filtres passe-bas, comme les 
integrales iterees de (4), laissent passer les basses frequences. Renvoyons a 
[4] pour lever ces contradictions apparentes. 

Remarque 5 Reprenons (2) et supprimons en (3) la condition tp < 1. Pour 
At — > et tp = N 2 = t^i, l a somme (3) ne tend plus presque surement vers 
0. Nos methodes sont done, par essence, non asymptotiques . 

Remarque 6 On ne pent quantifier, d'apres la proposition 3.2 et la remarque 
5, les performances d'un estimateur de type (4) dans le cadre habituel (cf. 
[15,21]), qui est probabiliste, asymptotique et hilbertien (fonctions de carres 
sommables) . Des criteres plus appropries seront proposes. 

3.2.2 Demodulation et erreurs en rafales 

En demodulation {cf. [17,21]), le parametre 9 a estimer, appele symbole, est 
associe, en general, a un signal, comme un sinus, cardinal ou non, solution 
d'une equation differentielle lineaire, a coefficients polynomiaux. Les erreurs en 
rafales, tres courantes, que les codes correcteurs d'erreurs doivent contrecarrer 
(cf. [17,21]), s'interpretent ainsi : l'estimation repose sur une hypothese fausse, 
a savoir le caractere centre du bruit. Comme la duree de transmission d'un 
symbole est « courte », on suppose que la moyenne m(t), a priori inconnue, 
est approchee par un polynome p{t) de degre limite, donne : p{t) est l'ombre 
de m(t). II semblerait que Ton puisse souvent supposer cette moyenne « a peu 
pres » constante pendant chaque duree de transmission, ce qui simplifie les 
calculs. 

Ecrivons n = n + p et considerons p comme une perturbation structured que 
Ton annihile [6,10] avec une puissance limitee, suffisante, de -|. On obtient une 
formule analogue a (4) en y remplagant n par no- La possible valeur du resultat 
suivant, de demonstration semblable a celle de la proposition 3.2, tient au fait 
que l'estimateur doit detecter un symbole parmi plusieurs, « eloignes » les 
uns des autres. 
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Proposition 3.3 Supposons que V ombre du signal ne soit pas annihilee par 
une puissance limitee de ~. II existe, alors, une estimee \0\ e (t) du symbole 
appartenant au halo de 6, pour toute largeur t n'appartenant pas au halo d'un 
zero du diviseur. 

Remerciements. L'auteur exprime sa reconnaissace a C. Lobry (Nice) et T. Sari 
(Mulhouse) pour d'utiles conversations. 
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